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ABSTRACT/RESUME 

Because  of  the  high  costs  associated  with  flight  and  experimental  measurements  in  wind 
tunnels  and  test  cells,  Computational  Fluid  Dynamics  (CFD)  is  becoming  a  useful  tool  for  predicting 
fluid  flow.  With  the  recent  performance  improvements  of  computers  and  numerical  codes,  a  much 
faster  prediction  of  fluid  flow,  shock  wave  position  and  aerodynamic  coefficients  calculation  is 
believed  to  be  possible  allowing  the  simulations  of  more  complex  flow. 

This  report  lies  in  the  continuity  of  exploring  phases  on  the  feasibility  to  use  Flow  Effectors  for 
Active  Flow  Control  for  the  Guidance  of  Missiles.  The  present  study  is  done  on  an  aerodynamic 
geometry  with  a  conical  nose  adapted  to  the  wind  tunnel  model.  Different  combinations  of  rectangular 
flow  effectors  place  on  the  nose  can  have  very  different  characteristics  in  steady  and  unsteady 
operation.  Each  one  of  these  items  has  a  direct  impact  on  the  general  performance  of  the  missile.  The 
missile  performances  can  also  be  affected  by  flight  parameters  such  as  Mach  number,  altitude  and 
angle  of  attack. 

The  aim  of  this  study  is  to  evaluate  the  performance  of  an  active  flow  control  in  the  form  of  a 
rectangular  flow  disturber  with  the  aim  of  controlling  a  missile.  Using  CFD,  combinations  of  different 
dimensions  of  flow  effectors  at  different  positions  on  the  surface  of  the  missile  were  evaluated  for 
their  impact  on  overall  system  performance  for  a  missile  flying  at  Mach  1.5  at  an  altitude  of  6  km  from 
zero  to  twenty  degrees  of  angle  of  attack. 


En  raison  des  couts  eleves  associes  aux  essais  en  vol  et  aux  tests  experimentaux  en  soufflerie, 
les  codes  numeriques  de  mecanique  des  fluides  sont  devenus  des  outils  tres  utiles  pour  predire  les 
ecoulements.  Avec  les  recents  developpements  dans  le  domaine  de  rinformatique  et  des  codes 
numeriques,  une  prediction  beaucoup  plus  rapide  de  I’ecoulement,  de  la  position,  de  la  propagation  des 
ondes  de  choc  et  de  revaluation  des  coefficients  aerodynamiques  est  rendue  possible,  permettant  ainsi 
la  simulation  de  problemes  plus  complexes. 

Ce  rapport  est  en  continuity  de  deux  phases  d’exploration  sur  la  possibility  d’utiliser  des 
perturbateurs  d’ecoulement  pour  le  Controle  Actif  d’Ecoulement  pour  le  Guidage  des  Missiles.  La 
presente  etude  est  executee  sur  une  geometric  aerodynamique  avec  un  nez  conique  adaptee  au  modele 
de  la  soufflerie.  Differentes  combinaisons  de  perturbateurs  rectangulaires  localises  sur  le  nez  peuvent 
avoir  differentes  caracteristiques  en  regime  permanent  et  transitoire.  Chacun  de  ces  items  a  un  impact 
direct  sur  la  performance  generate  du  missile.  Les  performances  du  missile  peuvent  aussi  etre 
affectees  par  les  parametres  de  vol  tels  que  le  nombre  de  Mach,  Paltitude  et  Tangle  d’attaque. 

Le  but  de  cette  etude  est  d’evaluer  le  rendement  d’un  controle  actif  d’ecoulement  sous  forme 
d’un  perturbateur  rectangulaire,  afin  de  controler  la  trajectoire  d’un  missile  ou  une  fusee  au  moyen 
d’un  Controle  Actif  d’Ecoulement.  En  utilisant  un  code  numerique  de  mecanique  des  fluides, 
differentes  dimensions  de  perturbateurs  avec  differentes  positions  sur  la  surface  du  missile  ont  ete 
evaluees  afin  de  verifier  leur  impact  sur  les  performances  globales  d’un  missile  volant  a  Mach  1.5  a 
une  altitude  de  6  km,  de  zero  a  vingt  degres  d’angle  d’attaque. 
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EXECUTIVE  SUMMARY 

The  DRDC  Valcartier  Flight  Mechanics  Group  is  responsible  for  research  in  aerodynamic  and  flight 
dynamic  behaviour  of  weapon  systems.  The  aeroballistic  facilities  are  being  used  to  obtain 
aerodynamic  coefficients  necessary  to  perform  trajectory  simulations  and  for  optimization  of  particular 
shapes  of  projectile  or  munitions  to  obtain  the  desired  aerodynamic  performance.  To  overcome  some 
of  the  physical  or  financial  constraints  imposed  by  these  experiments,  a  Computational  Fluid 
Dynamics  (CFD)  approach  can  be  used  to  complement  and  enhance  experimental  results. 

The  objective  of  the  contract  was  to  evaluate  the  performance  of  an  active  flow  control  in  the 
form  of  a  rectangular  flow  disturber  with  the  aim  of  controlling  a  missile  or  rocket.  Using  CFD, 
combinations  of  different  dimensions  of  flow  effectors  at  different  positions  on  the  surface  of  the 
missile  were  evaluated  for  their  impact  on  overall  system  performance  for  a  missile  flying  at  Mach  1.5 
at  wind  tunnel  conditions  from  zero  to  twenty  degrees  of  angle  of  attack. 

A  series  of  CFD  analysis  were  made  on  an  aerodynamic  geometry  with  some  Flow  Effectors. 
Three  different  dimensions  of  rectangular  Flow  Effectors  were  evaluated  for  their  side  force 
performance.  The  side  force  generated  by  the  Flow  Effectors  should  be  able  to  steer  a  missile.  For  the 
baseline  flow  effectors,  two  configurations  were  analysed.  A  configuration  with  three  3x  Flow 
effectors  were  calculated  and  two  configurations  with  five  times  the  baseline  Flow  Effectors  were 
calculated.  The  Flow  effectors  were  located  in  a  perpendicular  plane  with  the  angle  of  attack.  The 
experience  gained  during  the  previous  work  was  used  as  a  guideline  for  determining  the  models  and 
setting.  The  validation  was  done  in  a  previous  work  about  the  MEMS  feasibility  [DeChamplain  et  al. 
(2002)].  Finally,  recommendations  for  improvements  of  the  performance  were  stated. 

This  project  was  part  of  DRDC  technology  investment  found  (TIF)  project  “Supersonic  Missile 
Flight  Control  by  Manipulation  of  the  Flow  Structure  using  Micro-Actuated  Surfaces”.  It  resulted  in 
an  improved  capability  to  predict  numerically  the  aerodynamic  behaviour  of  rectangular  Flow 
Effectors  before  they  can  be  tested  in  DRDC  Valcartier  wind  tunnel. 


vii 


UNCLASSIFIED/  SANS  CLASSIFICATION 


1.0  INTRODUCTION 

Flow  effectors  (FE)  are  integrated  devices  located  in  strategic  places  on  the  surface  of  a 
missile.  These  systems  can  control  and  actuate  on  the  micro  scale,  and  function  individually  or  in 
arrays  to  generate  effects  on  the  macro  scale.  In  the  most  general  form,  MEMS  consist  of  mechanical 
microstructures,  micro  sensors,  micro  actuators  and  electronics  integrated  in  the  same  environment 
(i.e.  on  a  silicon  chip).  Miniaturization  of  mechanical  systems  promises  unique  opportunities  for  new 
directions  of  scientific  and  technological  progress. 

Missile  guidance  permits  an  increase  in  lethality  for  combat  superiority.  The  performance  of 
this  system  is  function  of  the  execution  quality  of  this  command.  The  majority  of  missiles  use 
aerodynamic  surfaces  to  accomplish  this  work.  The  control  of  the  thrust  vector  constitutes  another 
strategy  for  missile  control.  Nevertheless,  these  systems  require  a  structure  and  an  actuator.  This 
implies  number  of  systems  to  accomplish  the  task,  which  reduces  the  reliability  and  increases  the 
weight  of  the  mechanism.  Aerodynamic  surfaces  need  a  relatively  important  space  for  storage. 
Depending  on  the  desired  result,  one  might  wish  to  delay  or  accelerate  transition,  reduce  drag,  or 
increase  manoeuvrability.  Removing  conventional  controls  could  lead  to  significant  weight  reductions 
or  decrease  the  missile  signature. 

The  present  work  analyses  different  types  of  flow  effectors,  which  are  placed  in  a 
perpendicular  plane  with  the  angle  of  attack.  In  the  previous  work  with  MEMS,  [DeChamplain  et  al. 
(2001)]  and  [DeChamplain  et  al.  (2002)],  MEMS  devices  were  located  in  the  same  plane  of  the  angle 
of  attack.  Because  the  missile  when  at  an  angle  of  attack  hid  the  MEMS,  its  performance  was  found  to 
be  relatively  low.  A  general  objective  of  research  about  MEMS  is  to  investigate  the  magnitude  and 
type  of  surface  deformation  necessary  on  the  surface  of  a  missile  or  rocket  to  obtain  sufficient  control 
authority  for  effective  guidance  during  flight.  This  project  was  part  of  DRDC  technology  investment 
found  (TIF)  project  “Supersonic  Missile  Flight  Control  by  Manipulation  of  the  Flow  Structure  using 
Micro-Actuated  Surfaces”. 

The  aim  of  this  study  was  achieved  by  a  Computational  Fluid  Dynamics  (CFD)  analysis.  The 
commercial  CFD  code  Fluent™  was  used  to  accomplish  these  calculations.  One  should  note  that  the 
mesh  for  every  CFD  simulation  was  generated  with  the  commercial  code  Gambit™.  The  methodology 
used  included  many  steps.  The  experience  gained  during  the  previous  work  [DeChamplain  et  al. 
(2002)]  was  used  as  a  guideline.  In  the  previous  work,  a  good  validation  step  was  done  for  the 
geometry  without  MEMS. 

This  document  is  divided  in  many  parts.  At  first,  details  are  given  about  the  grid  and  the  model 
used  in  the  previous  work.  Furthermore  in  the  second  part,  a  description  of  the  flow  effectors,  their 
positions,  their  dimensions  and  the  different  angles  of  attack  that  were  used  are  described. 
Presentation  of  the  results  was  done  accordingly.  In  the  third  part,  an  analysis  of  the  results  is  given 
and  a  discussion  is  presented.  Finally,  a  conclusion  and  suitable  recommendations  are  given  in  the  last 
section. 
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2.0  METHODOLOGY 

When  a  Computational  Fluid  Dynamics  (CFD)  study  is  to  be  done,  a  complete  validation  study 
must  be  undertaken.  Basically,  this  means  that  CFD  results  must  be  verified  by  comparing  them  with 
experimental  data,  to  be  sure  that  the  CFD  model  can  predict  correctly  all  the  physical  phenomena’s. 
In  the  previous  study  on  MEMS  [DeChamplain  et  al.  (2002)]  a  good  validation  case  was  done.  The 
present  work  is  based  on  the  experience  gained  in  the  past.  The  information  about  the  various  models 
available  in  Fluent™  will  not  be  given  in  this  report.  Nonetheless,  it  could  be  useful  to  indicate  briefly 
some  parameters  that  were  selected  for  the  present  simulations.  Additionally,  a  validation  is  done  with 
a  mesh  that  is  refined  in  the  boundary  layer. 

2.1  Models 

Within  Fluent™  it  is  necessary  to  set  many  parameters  for  calculation  purposes.  All  the 
possibilities  will  not  be  described  here.  Only  parameters  chosen  for  the  present  analysis  will  be 
mentioned.  Readers  who  would  like  more  details  about  those  parameters  should  consult 
[DeChamplain  et  al.  (2001)]  or  Fluent™  literature. 

In  the  present  simulations,  Navier-Stokes  equations  are  resolved  with  the  coupled  form  of  the 
code  where  the  energy  equation  is  solved  simultaneously  with  the  fluid  equations.  Those  equations 
are  solved  with  an  explicit  pseudo-time  formulation.  This  formulation  calculates  more  rapidly  and 
need  less  computer  memory.  The  k-e  model  improved  with  the  Realizable  model  is  adopted  to  take 
into  account  turbulence  and  vortex  formation.  Finally,  the  effects  of  boundary  layer  are  evaluated  with 
the  Enhanced-Wall  treatment  model.  This  model  proposes  two  options  to  improve  the  calculation.  The 
pressure  one  is  active,  but  since  the  thermal  heat  transfer  is  not  calculated,  the  thermal  option  is  not 
activated.  The  working  fluid  is  air.  The  density  is  calculated  with  the  ideal-gas  law  and  the  viscosity, 
with  the  Sutherland  three  coefficients  model.  All  the  other  properties  are  as  the  default  settings. 

2.2  Grids 

In  the  present  study  it  is  not  possible  to  split  the  calculation  domain  in  two  parts  like  it  was 
done  in  the  previous  work  [DeChamplain  et  al.  (2002)].  Since  the  Flow  effectors  (FE)  are  located  in  a 
perpendicular  plane  with  the  angle  of  attack,  the  symmetry  plane  that  existed  in  the  previous  work 
disappeared  in  this  case.  The  complete  domain  was  calculated  in  the  present  3D  calculation.  As  the 
present  geometry  had  a  small  rounded  nose,  it  was  easier  to  build  mesh  upstream  of  the  geometry.  A 
small  cylinder  was  built  starting  from  the  nose  to  the  upstream  part  of  the  domain.  The  diameter  of  the 
cylinder  was  around  0,003  inches  compared  to  the  0,005  inches  round  tip  of  the  nose  reducing  the 
mesh  skewness  on  the  nose  tip. 

In  accordance  with  the  previous  MEMS  study,  grids  generated  were  structured.  It  has  been 
demonstrated  that  an  unstructured  mesh  shows  instabilities  and  cannot  give  reasonable  solutions. 
However,  some  differences  exist  in  the  mesh  employed  presently.  All  of  them  are  related  to 
rectangular  flow  effectors  that  needed  some  refinement  in  the  grid  in  their  vicinity.  Since  the 
Enhanced-Wall  treatment  was  used,  the  mesh  had  to  have  a  Y+  parameter  around  1  over  all  the  missile 
wall.  The  supersonic  shock  wave  presence  makes  it  difficult  to  establish  the  Y+  with  the  same  value 
over  the  entire  missile.  An  effort  was  made  to  have  a  more  favourable  Y+  on  the  Flow  effectors  (FE) 
device  and  region.  In  the  results  tables,  the  Y+  is  calculated  over  the  entire  missile  and  around  the  FE 
region.  It  is  to  be  considered  as  a  mean  Y+.  It  may  be  important  to  remember  that  the  Y+  parameter  is 
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calculated  like  a  Reynolds  number  at  the  first  cell  near  the  wall.  Equation  (1)  is  used  to  calculate  the 
Y+  value.  This  value  is  normally  calculated  by  Fluent™  code.  This  equation  includes  density  pw, 
velocity  Vw,  height  of  the  first  cells  Lw  and  viscosity  /Av-  All  those  values  are  calculated  at  the  first 
cell  near  the  wall. 

The  Flow  effectors  (FE)  analysed  in  the  present  study  had  the  shape  of  a  rectangular  effectors. 
Three  FE  were  analyzed,  the  baseline  effectors,  the  three  time  baseline  effectors  and  the  five  time 
baseline  effectors.  The  relative  position  and  number  of  FE  were  the  parameters  of  interest  for  each  FE 
dimension.  The  calculation  included  1,  3  or  5  FE.  All  FE  were  located  at  25,32  mm  (0,997  inches) 
from  the  nose  tip.  The  space  between  the  FE  was  defined  by  an  angle  around  the  axis  of  the  geometry. 
The  angle  between  the  FE  was  15  or  30  degrees  and  for  the  3  FE  geometries,  the  angle  was  45 
degrees.  The  angles  of  attack  for  the  study  were  0,  10,  15  and  20  degrees. 

A  view  of  the  mesh  in  the  far  field  is  presented  in  figure  A-l.  Details  of  the  refinement 
necessary  around  the  FE  are  presented  in  figures  A-2  and  A-3  with  the  baseline  FE.  To  simplify  the 
meshing  around  the  FE,  the  domain  was  split  into  volumes.  Structured  volumes  are  presented  in 
figure  A-4  for  the  baseline  FE,  in  figure  A-5  for  the  3x  FE  and  in  figure  A-6  for  the  5x  FE.  Finally, 
figures  A-l  and  A-8  give  an  example  of  the  mesh  around  the  3x  FE.  Figures  A-l  to  A-8  are  presented 
in  appendix  A. 


Y+  _ 


[1] 


2.3  Boundary  Conditions 

The  domain  is  surrounded  by  a  pressure  far-field  boundary  condition  in  front  and  around  the 
calculation  domain.  The  outlet  of  the  domain  determined  by  a  pressure  outlet  boundary  condition  to 
reduce  the  boundary  condition  effect  and  to  be  sure  that  no  shock  wave  reflection  can  exist.  The 
missile  is  flying  according  to  new  conditions  compared  to  the  previous  MEMS  study.  The  missile  has 
a  design  Mach  number  of  1,5.  The  static  temperature  is  206  K  and  the  static  pressure  is  26  000  Pa. 
Furthermore,  three  angles  of  attack  (AO A)  were  selected  for  calculations:  10  and  15  degrees.  Angle  of 
attack  of  0  and  20  degrees  were  kept  as  an  option  if  time  permitted. 

2.4  Aerodynamic  Coefficients 

The  results  presented  will  the  aerodynamic  coefficients.  The  aerodynamic  coefficients  are 
calculated  with  the  missile  coordinate  axis.  The  origin  of  the  coordinate  axis  is  located  at  the  nose  tip. 
The  aerodynamic  moment  is  calculated  at  the  base  of  the  missile.  Figure  1  presents  the  coordinate 
system  for  the  aerodynamic  coefficients  and  the  moment  center  for  the  moment  calculation.  The  base 
drag  is  not  calculated  in  the  x  axis  aerodynamic  coefficient.  Table  1  presents  the  reference  values  used 
in  the  calculations. 
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Figure  1:  Coordinate  system  for  the  calculation  of  the  aerodynamic  coefficients 

and  moment. 


Area  (  A ),  m2 

706,9  E-6 

Temperature  (  T),  K 

206 

Density  ( p),  kg/m3 

0,439  7104 

Velocity  (  V),  m/s 

431,423 

Enthalpy  (  h ),  J/kg 

300  387,8 

Viscosity  (  jj  ),  kg/m  s 

1,789  E-5 

Length  (  L ),  m 

0,03 

Ratio  of 
specific  heats 

1,4 

Pressure  (  P ),  Pa 

26  000 

Table  1  Reference  values  for  the  calculation  of  the  aerodynamic 

COEFFICIENTS 
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3.0  CALCULATIONS 

To  evaluate  the  effect  of  the  FE  on  the  flow  and  the  accuracy  of  the  calculation,  the  body  alone 
geometry  was  simulated  at  every  angle  of  attack.  The  results  were  compared  with  Missile  Datcom 
97™  results  to  have  an  idea  of  the  length  scale  of  the  solution.  The  Missile  Datcom  97™  calculation 
had  a  difference  of  around  25%  with  the  wind  tunnel  results.  The  body  alone  geometry  is  shown  in 
figure  2.  As  presented  in  figure  1,  the  x  axis  is  on  the  centre  axis  of  the  geometry.  The  angle  of  attack 
is  in  the  X-Y  plane  and  the  FE  devices  are  located  on  the  positive  Z  side  of  the  geometry.  The 
diameter  of  the  geometry  is  30  mm.  The  conical  part  of  the  geometry  is  90  mm  in  length  and  the 
cylindrical  part  is  300  mm  in  length.  The  total  length  of  the  geometry  is  390  mm  (15,353  inches). 


Figure  2:  Body  alone  geometry  configuration  with  the  3D  coordinate  axis 


3.1  Flow  effectors  dimensions 

The  Flows  effectors  (FE)  dimensions  are  simpler  compared  to  the  microbubbles  of 
[DeChamplain  et  al.  (2002)].  The  length  of  the  FE  was  defined  in  the  x  axis  of  the  aerodynamic 
geometry.  The  width  of  all  FE  is  0,78  mm  (0,031  inches).  The  baseline  FE  was  approximately  1  mm 
(0,0394  in)  high  and  2,67  mm  long  (0,105  in).  The  three  time  baseline  (3x)  FE  were  near  3  mm  by  8 
mm  (0,1 18  in  by  0,315  in)  and  the  five  time  baseline  (5x)  FE  were  around  8  mm  by  13,3  mm  or  (0,197 
in  by  0,515  in).  The  exact  dimensions  of  all  FE  are  presented  in  figure  3.  The  FE  dimensions  are  all 
defined  in  inches  on  the  figures. 
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Figure  6:  Position  of  the  5x  FE  on  the  aerodynamic  geometry 

3.3  Results  without  Flow  effectors 

To  begin  with,  a  calculation  of  the  aerodynamic  geometry  was  done  to  define  a  mesh  having  a 
Y+  parameter  near  one.  The  first  mesh  gave  a  Y+  value  near  0,03  in  the  Flow  effectors  (FE)  zone. 
With  a  too  small  value  of  Y+,  a  second  calculation  presented  a  mean  value  of  Y+  of  0,45  in  the  FE 
zone.  Since  no  experimental  data  were  available  to  verify  the  exactitude  of  the  calculation  when  the 
study  begun,  it  was  decided  to  build  another  mesh  to  having  a  Y+  value  more  near  of  the  ideal  value  of 
one.  The  last  mesh  had  an  Y+  mean  value  in  the  FE  zone  of  0,88.  Along  the  strait  part  of  the  body, 
the  flow  is  less  aggressive  on  the  geometry  wall,  and  the  calculated  Y+  value  is  lower.  The  Y+  over  all 
the  aerodynamic  geometry  varied  between  0,6  and  0,67  with  an  angle  of  attack.  The  height  of  the  first 
cells  was  at  1,5  x  10'4  inches  to  obtain  the  actual  Y+  value.  The  calculation  for  all  the  study  included 
two  values  of  Y+.  The  first  one  was  the  mean  value  in  the  FE  zone  and  the  FE  itself.  The  second  one 
was  the  mean  value  over  all  the  aerodynamic  body  and  it  included  the  FE  itself  for  the  calculation  with 
the  FE.  Table  2  include  the  condition  of  the  calculation  without  FE.  Table  3  include  the  aerodynamic 
coefficients  results  for  the  calculation  without  the  FE.  The  results  for  Missiles  Datcom™  calculation 
are  also  included  in  the  table.  Without  FE,  CY  should  be  zero;  this  value  in  the  table  gives  an  idea  of 
the  precision  of  the  calculation.  FE  side  force  under  the  present  value  may  represent  a  negligible  effect 
of  the  FE. 
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Name 

Mesh 

AOA 

Nb  cells 

Y+  (FE)  i  Y+  (All) 

Without  FE 

Mis-CC_AOAOO 

893276 

0,87 

0,60 

0 

Mis-CC_AOA10 

893276 

0,88 

0,64 

10 

Mis-CC_AOA15 

893276 

0,88 

0,66 

15 

Mis-CC_AOA20 

893276 

0,88 

0,67 

20 

Missile  Datcom 

Datcom  AOAOO 

... 

0 

Datcom  AOA1 0 

10 

Datcom  AOA1 5 

--- 

... 

15 

Datcom  AOA20 

— 

... 

20 

Table  2  Numerical  configuration  without  FE 


Name 

Forces 

Moments 

Cx 

Cn 

Cy 

c, 

C„ 

C  m 

Without  FE 

Mis-CC_  AOAOO 

0,2399 

0,0000 

0,0001 

0,0000 

0,0004 

0,0001 

Mis-CC_AOA10 

0,2429 

0,6601 

0,0004 

0,0000 

0,0016 

-5,7154 

Mis-CC_AOA15 

0,2309 

1 ,3234 

0,0087 

-0,0001 

0,0306 

-10,6938 

Mis-CC_AOA20 

0,2084 

2,5469 

0,0512 

-0,0001 

0,2422 

-19,4643 

Datcom 

Datcom  AOAOO 

0,2310 

0,0000 

0,0000 

0,0000 

0,0000 

0,0000 

Datcom  AOA10 

0,2027 

0,9780 

0,0000 

0,0000 

0,0000 

4,8660 

Datcom  AOA1 5 

0,1691 

1,5760 

0,0000 

0,0000 

0,0000 

8,1840 

Datcom  AOA20 

0,1238 

3,2670 

0,0000 

0,0000 

0,0000 

19,4970 

Table  3  Numerical  coefficients  results  without  FE 


AOA 

Aerodynamic  coefficient 

Cx 

Cn 

C  m 

0 

3,8% 

... 

... 

10 

19,8% 

-32,5% 

17,5% 

15 

36,5% 

-16,0% 

30,7% 

20 

68,3% 

-22,0% 

-0,2% 

Table  4  Error  between  numerical  results  Missile  Datcom  97™  results 


I 

I 
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3.4  Results  with  Flow  effectors 

Since  the  calculations  without  Flow  effectors  (FE)  showed  acceptable  results,  it  was  now 
possible  to  calculate  the  performance  of  the  FE  with  good  confidence.  The  results  are  included  in  two 
tables.  Table  5  presents  the  calculation  conditions  and  settings.  This  table  includes  the  number  of 
cells  in  the  mesh  and  two  mean  Y+  values.  The  table  includes  FE  numbers  and  their  positions. 
Finally,  the  table  presents  the  angle  of  attack  of  calculation.  Table  6  includes  the  performance  results 
as  aerodynamic  coefficients  and  moments. 


Name 

Mesh 

Flow  effectors 

AOA 

Nb  cells 

Y+(FE) 

Y+  (All) 

Number 

_ type _ 

Baseline  (FE) 

MisC-A1_AOA10 

1129256 

0,90 

0,65 

1 

1 

10 

MisC-A2_AOA1 0 

1125750 

0,86 

0,64 

3 

2x  15° 

10 

MisC-A3_AOAOO 

1125750 

0,85 

0,60 

3 

2x  30° 

0 

MisC-A3_AOA1 0 

1125750 

0,88 

0,64 

3 

2x  30° 

10  I 

MisC-A3_AOA1 5 

1125750 

0,88 

0,66 

3 

2x  30° 

15 

MisC-A3_AOA20 

1125750 

0,88 

0,66 

3 

2x  30° 

20 

MisC-A4_AOA1 0 

1104946 

0,88 

0,64 

3 

2x  45° 

10 

MisC-A5_AOAOO 

1293010 

0,83 

0,60 

5 

4x  15° 

0 

MisC-A5_AOA1 0 

1293010 

0,85 

0,64 

5 

4x  15° 

10 

MisC-A5_AOA1 5 

1293010 

0,85 

0,66 

5 

4x  15° 

15 

MisC-A5_AOA20 

1293010 

0,85 

0,66 

5 

4x  15° 

20 

[Three  time  baseline  (3x  FE) 

MisC-B3_AOAOO 

1142698 

0,85 

0,60 

3 

2x  30° 

0 

MisC-B3_AOA10 

1142698 

0,88 

0,65 

3 

2x  30° 

10 

MisC-B3_AOA15 

1142698 

0,89 

0,66 

3 

2x  30° 

15 

MisC-B3_AOA20 

1142698 

0,89 

0,67 

3 

2x  30° 

20 

Five  time  baseline  (5x  FE) 

MisC-C1_AOAOO 

1300894 

0,93 

0,61 

1 

0 

MisC-C1_AOA10 

1300894 

0,94 

0,65 

1 

1 

10 

MisC-C1_AOA15 

1300894 

0,95 

0,66 

1 

1 

15 

MisC-C1_AOA20 

1300894 

0,94 

0,67 

1 

1 

20 

M  isC-C4_AOAOO 

1385693 

0,88 

0,60 

3 

2x  45° 

0 

MisC-C4_AOA10 

1385693 

0,89 

0,64 

3 

2x  45° 

10  ■ 

MisC-C4_AOA15 

1385693 

0,89 

0,66 

3 

2x  45° 

15 

M  isC-C4_AOA20 

1385693 

0,90 

0,67 

3 

2x  45° 

20 

Table  5  Numerical  configuration  with  FE 
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Name 

Results 

Cx 

CN 

Cy 

c, 

c„ 

C  m 

Baseline  (FE) 

MisC-A1_AOA10 

0,2450 

0,6602 

-0,0397 

-0,0002 

-0,1559 

-5,7326 

MisC-A2_AOA10 

0,2463 

0,6598 

-0,0353 

-0,0003 

-0,1413 

-5,7322 

MisC-A3_AOAOO 

0,2428 

0,0000 

-0,0007 

0,0000 

-0,0052 

-0,0001 

MisC-A3_AOA10 

0,2469 

0,6617 

-0,0777 

-0,0003 

-0,3137 

-5,7475 

MisC-A3_AOA1 5 

0,2352 

1,3415 

-0,3439 

-0,0006 

-1,4170 

-10,7895 

MisC-A3_AOA20 

0,2132 

2,5928 

-0,4648 

-0,0017 

-2,5391 

-19,6688 

MisC-A4_AOA1 0 

0,2472 

0,6634 

-0,0998 

-0,0003 

-0,4025 

-5,7509 

MisC-A5_AOAOO 

0,2443 

0,0000 

-0,0012 

0,0000 

-0,0097 

0,0000 

MisC-A5_AOA1 0 

0,2485 

0,6596 

-0,0362 

-0,0004 

-0,1509 

-5,7308 

MisC-A5_AOA1 5 

0,2365 

1 ,3250 

-0,1586 

-0,0005 

-0,6230 

-10,7189 

MisC-A5_AOA20 

0,2149 

2,5637 

-0,3107 

-0,0012 

-1,6377 

-19,5481 

Three  time  baseline  (3x  FE) 

MisC-B3_AOAOO 

0,2528 

0,0000 

-0,0021 

0,0000 

-0,0141 

-0,0001 

MisC-B3_AOA10 

0,2584 

0,6919 

-0,2290 

-0,0045 

-1,0811 

-6,0329 

MisC-B3_AOA1 5 

0,2469 

1,3900 

-0,5686 

-0,0077 

-2,9107 

-11,2638 

M  isC-B3_AO  A20 

0,2268 

2,5983 

-0,4338 

-0,0103 

-3,1138 

-20,0635 

Five  time  baseline 

(5x  FE) _ 

MisC-C1_AOA00 

0,2490 

0,0000 

-0,0007 

0,0000 

-0,0003 

-0,0001 

MisC-C1_AOA10 

0,2529 

0,7118 

-0,2566 

-0,0107 

-1 ,3651 

-6,2976 

MisC-C1_AOA15 

0,2430 

1,3947 

-0,4430 

-0,0177 

-2,9600 

-11,5818 

MisC-C1_AOA20 

0,2235 

2,6862 

0,0037 

-0,0234 

-2,0535 

-20,7208 

MisC-C4_AOAOO 

0,2638 

0,0000 

-0,0030 

0,0000 

-0,0199 

-0,0002 

MisC-C4_AOA10 

0,2673 

0,7385 

-0,3901 

-0,0144 

-2,0864 

-6,4674 

MisC-C4_AOA15 

0,2565 

1,4580 

-0,6996 

-0,0237 

-4,4060 

-1 1 ,9883 

MisC-C4_AOA20 

0,2390 

2,7342 

-0,2779 

-0,0338 

-4,0068 

-21 ,2730 

Table  6  Results  of  the  performance  of  the  FE  system 

It  is  important  to  consider  that  a  CY  value,  which  is  the  side  force  of  the  FE  device,  that  has  a 
value  below  the  calculations  without  FE  has  a  negligible  effect.  Table  6  shows  a  negligible  effect  of 
the  FE  with  an  angle  of  attack  of  0  degrees.  With  an  angle  of  attack  of  10  degrees,  the  side  force  was 
around  one  tenth  of  the  drag  for  baseline  FE.  The  side  force  was  bigger  than  the  drag  for  5x  FE. 
Normally  the  performance  of  the  FE  increases  when  the  angle  of  attack  increases.  But,  with  an  angle 
of  attack  of  20  degrees,  the  boundary  layer  was  separating  over  the  FE,  reducing  the  side  force  of  the 
devices.  For  the  calculation  with  a  5x  FE  at  20  degrees  of  angle  of  attack,  the  side  force  effect  was 
changing  of  side.  Generally,  side  force  applied  on  the  missile  by  a  FE  stirred  the  missile  in  the 
opposite  side  of  its  original  position. 

Fluent™  solutions  are  presented  on  some  contours  in  appendix  B.  Figures  B-l  and  B-2  present 
a  general  flow  field  with  pressure  contours  in  plane  X-Y,  with  an  angle  of  attack,  and  plane  X-Z,  with 
the  FE.  These  two  figures  were  from  the  3x  FE  with  an  angle  of  attack  of  10  degrees.  The  shock 
waves  from  the  geometry  with  an  angle  of  attack  are  easy  to  see  in  figure  B-l.  Effects  of  FE  are 
presented  in  figure  B-2.  Figure  B-3  shows  the  Mach  number  around  three  baselines  FE  and  figure  B-4 
presents  the  Mach  number  around  five  baselines  FE.  These  two  figures  help  understanding  why  side 
force  performance  is  lower  when  FE  having  an  angle  of  15  degrees  between  them.  Figure  B-4  shows 
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that  the  Mach  number  stays  small  between  the  FE.  When  the  space  between  FE  is  increased,  as  in 
figure  B-3,  the  flow  can  pass  through  all  FE  without  being  perturbed.  A  special  effect  was  the 
performance  lost  that  happened  with  the  3x  FE  and  5x  FE  at  an  angle  of  attack  of  20  degrees.  To 
understand  this  phenomenon,  Mach  number,  static  pressure  field  and  dynamic  pressure  field  were 
compared.  Figures  B-5  to  B-10  compare  solutions  with  3x  FE  with  an  angle  of  attack  of  15  degrees 
and  20  degrees.  Mach  number  and  static  pressure  fields  were  presented  around  the  FE.  Dynamic 
pressure  field  was  captured  at  the  half-length  of  the  missile  to  see  if  the  generated  vortex  was 
perturbed.  Figures  with  Mach  number  contours  present  a  high  velocity  over  the  upper  FE  with  an 
angle  of  attack  of  15  degrees.  With  an  angle  of  attack  of  20  degrees,  a  boundary  layer  separation 
occurred  over  the  upper  FE.  A  small  vortex  created  a  low  velocity  zone  over  the  lower  FE  with  an 
angle  of  attack  of  15  and  20  degrees.  These  two  phenomena’s  can  explain  the  low  side  force  ay  20 
degrees  AOA.  Static  pressure  around  the  two  solutions  had  small  differences.  A  pressure  difference 
was  calculated  under  the  missile  because  of  the  effect  of  the  variation  of  angle  of  attack.  The 
boundary  layer  separation  had  the  effect  that  the  pressure  stayed  higher  with  an  angle  of  attack  of  20 
degrees.  The  figure  with  the  dynamic  pressure  contours  shows  that  the  3x  FE  modified  the  vortex 
shape  created  by  the  missile.  The  difference  between  the  two  calculations  were  small  and  it  did  not 
appeared  to  explain  the  lost  of  side  force  when  the  angle  of  attack  was  increased  to  20  degrees. 
Finally,  the  flow  field  stayed  very  symmetric  around  the  FE  when  the  angle  of  attack  was  zero.  Figure 
B-l  1  was  calculated  with  three  5x  FE.  The  small  performance  of  the  FE  without  an  angle  of  attack 
could  be  explained  by  the  fact  that  the  FE  were  hidden  by  the  shock  wave  and  by  the  symmetric  flow 
around  the  FE. 

Three  graphs  were  built  for  side  force  performance  of  the  FE  with  an  angle  of  attack.  In  those 
graphs,  a  side  force  is  calculated  from  15  degrees  AOA  on  the  missile  with  no  flow  effectors.  This 
could  be  explained  from  the  fact  that  the  mesh  is  not  equally  distributed  around  the  body  of  the  missile 
and  that  small  mesh  distortions  could  induce  asymmetrical  vortices,  which  could  induce  a  side  force 
on  the  missile  body.  Figure  C-l  presents  the  performance  of  a  baseline  FE.  Figure  C-2  presents  the 
side  force  of  3x  FE.  .  For  the  3x  FE,  we  can  see  that  the  side  force  drops  from  15  to  20  degrees  AOA. 
This  is  explained  by  the  fact  that  the  upper  flow  effector  is  in  a  stall  on  its  upper  surface  as  seen  in 
figures  B-5  to  B-8  in  appendix  B.  Figure  C-3  presents  the  performance  of  the  5x  FE.  With  one  5x 
FE,  the  low  side  force  performance  can  be  explained  by  a  boundary  layer  separation  near  the  body  of 
the  missile.  For  the  configuration  with  three  FE,  some  vortex  were  generate  at  the  tip  of  the  FE.  The 
vortex  generated  a  low  velocity  region  in  the  middle  of  the  FE  with  an  angle  of  attack  of  15  and  20 
degrees.  The  low  velocity  zone  create  by  the  vortex  took  a  more  important  area  with  an  angle  of 
attack  of  20  degrees,  which  explain  the  lower  performance.  On  those  graphs,  side  force  is  presented  as 
a  pushing  effect  on  the  missile  to  obtain  positive  value.  Positive  values  are  easier  to  compare  the  side 
force  performance  of  the  FE. 

3.5  Parametric  study  with  Flow  effectors 

With  baseline  Flow  effectors  (FE),  some  additional  configurations  were  calculated.  Those 
calculations  permitted  to  understand  the  effect  of  the  angle  between  FE  and  the  number  of  FE.  For 
those  two  parametric  studies,  the  angle  of  attack  was  10  degrees.  In  the  first  study,  effect  of  the  angle 
between  FE  was  analysed.  The  results  are  presented  on  graphs  C-4  in  appendix  C.  Results  show  that 
when  the  angle  between  FE  was  increased,  performance  increased.  But  results  also  show  that  at  an 
angle  of  15  degrees  between  the  FE  was  not  enough  to  have  an  effect  on  side  force  performance.  Side 
force  approximately  doubled  from  one  FE  when  the  FE  had  an  angle  of  30  degrees  between  them.  The 
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side  force  was  2,5  times  bigger  with  FE  having  an  angle  of  45  degrees  between  them,  than  with  only 
one  FE. 

The  second  parametric  study  concerned  the  number  of  FE.  The  FE  that  had  an  angle  of  15 
degrees  between  them  made  it  easier  the  parametric  study  with  this  angle.  The  study  compares  the 
side  force  without  FE,  with  one,  three  and  five  FE.  Results  are  presented  on  graphs  C-5  in  appendix 
C.  The  graph  presents  a  small  variation  of  the  side  force  with  the  number  of  FE. 
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4.0  CONCLUSION 

This  study  on  flow  effectors  (FE)  aimed  at  defining  the  influence  of  the  number  and  the 
position  of  the  FE  on  a  missile  body  and  the  effect  of  the  angle  of  attack.  Furthermore,  three 
dimensions  of  rectangular  FE  were  analysed.  It  was  desired  to  know  if  a  given  configuration  for  FE 
could  have  a  greater  impact  than  others.  Since  the  validation  was  completed  in  previous  work,  the 
model  selection  was  defined  with  the  experience  gained.  From  the  results  obtained  with  numerical 
simulations,  some  remarks  can  be  given.  First  of  all,  the  FE  can  achieve  a  better  side  force  compared 
to  the  ones  calculated  in  the  previous  work.  The  performance  is  very  sensitive  to  the  angle  of  attack. 
When  the  angle  of  attack  is  zero,  the  side  force  is  near  zero.  The  side  force  coefficient  of  the  FE  can 
be  as  two  times  the  drag  for  baseline  FE  for  an  angle  of  attack  of  20  degrees.  The  side  force,  for  the 
3x  FE,  is  more  than  two  times  the  drag  at  an  angle  of  attack  of  15  degrees,  but  the  performance 
reduces  when  the  angle  of  attack  is  increased  to  20  degrees.  The  5x  FE  also  permits  an  important  side 
force  at  an  angle  of  attack  of  15  degrees,  but  an  important  performance  lost  happen  when  the  when  the 
angle  of  attack  is  20  degrees.  For  one  5x  FE  with  an  angle  of  attack  of  20  degrees,  the  side  force 
becomes  near  zero  and  the  force  is  switching  side.  The  FE  creates  an  important  force  that  can  steer  the 
missile  in  the  opposite  direction  of  its  original  position.  In  a  future  work,  it  will  be  interesting  to  have 
a  more  complete  parametric  study  about  the  number  of  FE  having  a  different  angle  between  them. 
Another  parametric  study  can  be  done  to  identify  the  effect  of  the  dimensions  of  the  FE  and  the 
optimal  dimensions.  The  present  study  presents  a  good  performance  of  the  FE  when  they  are  placed  in 
a  perpendicular  plane  with  an  angle  of  attack.  The  previous  work  had  resulted  in  a  poor  performance 
with  the  FE  in  the  same  plane  with  the  angle  of  attack.  Another  question  remains,  what  would  be  the 
performance  of  the  FE  with  an  intermediate  plane  with  the  angle  of  attack  ? 
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APPENDIX  A 


MESH  FIGURES  WITH 
RECTANGULAR  FLOW  EFFECTORS 
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Figure  A-l  View  of  the  mesh  in  the  far  Field  and  boundary  conditions  position 


Figure  A-2  Zoomed  view  of  the  mesh  around  the  baseline  FE 
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Figure  A-3  Boundary  refinement  around  the  baseline  FE 
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Figure  A-4  Block  construction  around  the  baseline  FE 
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Figure  A-5  Block  construction  around  the  3x  FE 


Figure  A -6  Block  construction  around  the  5x  FE 
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Figure  A-7  Mesh  cut  around  the  3x  FE  in  the  X-Z  plane 


Figure  A-8  Mesh  cut  around  the  3x  FE  in  the  Y-Z  plane 
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APPENDIX  B 


FLUENT™  SOLUTION  FIGURES  WITH 
RECTANGULAR  FLOW  EFFECTORS 
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Figure  B-l  Pressure  field  in  the  plane  of  the  angle  of  attack  with  three  5x  FE  at  10°  AOA 
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Figure  B-2  Pressure  field  in  the  plane  of  the  three  5x  FE  at  10°  AOA 
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FLUENT  6.1  (3d,  coupled  exp,  rke) 


Figure  B-3  Mach  number  around  3  baseline  FE  with  30  degrees  between  them  with  20°  AOA 
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FLUENT  6.1  (3d,  coupled  exp,  rke) 


Figure  B-4  Mach  number  around  5  baseline  FE  with  15  degrees  between  them  with  20°  AOA 


22 


UNCLASSIFIED/  SANS  CLASSIFICATION 


2  300+00 
2  21o+00 
2  140+001 
2  070+00 
2  00o+ 

1  930+ 

1  860+ ... 
1  79e+00| 
1  7  2e+“ 

1  660+ 

1  590+ 

1  52o+ 

1  45e+ 

1  38e+ 

1  310+ 

1  24e+ 

I  lOet 

1  03e+  _ 
Be-0 
7e-01 
“  e— 0 1 

__e-01 
900-01 
>210-01 

5  520-01 
4  830-01 
4 140-01 
3450-01 

2  760-01 
2  07< 

6  90e-02 
0  00e+00 


Contours  of  Mach  Number 


Mar  25,  2004 
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Figure  B-5  Mach  number  around  3  3x  FE  with  30  degrees  between  them  with  15°  AOA 
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FLUENT  6.1  (3d,  coupled  exp,  rke) 

Figure  B-6  Mach  number  around  3  3x  FE  with  30  degrees  between  them  with  20°  AOA 
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Figure  B-7  Pressure  around  3  3x  FE  with  30  degrees  between  them  with  15°  AOA 
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FLUENT  6.1  (3d,  coupled  exp,  rke) 


Figure  B-8  Pressure  around  3  3x  FE  with  30  degrees  between  them  with  20°  AOA 
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Mar  25,  2004 
FLUENT  6.1  (3d,  coupled  exp,  rke) 


Figure  B-9  Dynamic  pressure  in  the  middle  of  the  missile  (0,2  m)  with  15°  AOA  in  the  calculation 
with  3  3x  FE  having  an  angle  of  30  degrees  between  them 


Contours  of  Dynamic  Pressure  (pascal)  Mar  25,  2004 

FLUENT  6.1  (3d,  coupled  exp,  rke) 


Figure  B-10  Dynamic  pressure  in  the  middle  of  the  missile  (0,2  m)  with  20°  AOA  in  the  calcu 
with  3  3x  FE  having  an  angle  of  30  degrees  between  them 
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Contours  ot  Mach  Number  Mar  25,  2004 

FLUENT  6.1  (3d,  coupled  exp,  rke) 


Figure  B-l  1  Mach  number  around  a  5x  FE  with  0°  AOA 
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APPENDIX  C 


GRAPH  OF  THE  PERFORMANCE  OF 
RECTANGULAR  FLOW  EFFECTORS 
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igure  C-l  Graph  of  the  performance  of  the  baseline  FE  vs  the  angle  of  attack 


Figure  C-2  Graph  of  the  performance  of  the  3x  FE  vs  the  angle  of  attack 
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Angle  of  attack  in  degrees 

Figure  C-3  Graph  of  the  performance  of  the  5x  FE  vs  the  angle  of  attack 


Figure  C-4  Graph  of  the  performance  of  the  baseline  FE  vs  the  angle  between  the  FE  for  an  angle  of 
attack  is  10  degrees.  The  number  of  FE  is  limited  to  three.  With  an  angle  between  the  devices  of  zero 

degrees,  there  is  only  one  FE. 
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Figure  C-5  Graph  of  the  performance  of  the  baseline  FE  vs  the  number  of  FE  having  an  angle  of  15 
degrees  between  them  and  at  an  angle  of  attack  of  10  degrees 
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Because  of  the  high  costs  associated  with  flight  and  experimental  measurements  in  wind  tunnels  and  test  cells,  Computational 
Fluid  Dynamics  (CFD)  is  becoming  a  useful  tool  for  predicting  fluid  flow.  With  the  recent  performance  improvements  of 
computers  and  numerical  codes,  a  much  faster  prediction  of  fluid  flow,  shock  wave  position  and  aerodynamic  coefficients 
calculation  is  believed  to  be  possible  allowing  the  simulations  of  more  comolex  flow. 

This  report  lies  in  the  continuity  of  exploring  phases  on  the  feasibility  to  use  Flow  Effectors  for  Active  Flow  Control  for 
the  Guidance  of  Missiles  or  Rockets.  The  present  study  is  done  on  a  aerodynam  c  geometry  with  a  conical  nose  adapted  to  the 
wind  tunnel  model.  Different  combinations  of  rectangular  flow  effectors  place  on  the  nose  can  have  very  different  characteristics 
in  steady  and  unsteady  operation.  Each  one  of  these  items  has  a  direct  impact  on  the  general  performance  of  the  missile.  The 
missile  performances  can  also  be  affected  by  flight  parameters  such  as  Mach  number,  altitude  and  angle  of  attack. 

The  aim  of  this  study  is  to  evaluate  the  performance  of  an  active  flow  control  in  the  form  of  a  rectangular  flow  disturber  with  the 
aim  of  controlling  a  missile  or  rocket.  Using  CFD,  combinations  of  different  dimensions  of  flow  effectors  at  different  positions  on 
the  surface  of  the  missile  were  evaluated  for  their  impact  on  overall  system  performance  for  a  missile  flying  at  Mach  1.5  at  an 
altitude  of  6  km  from  zero  to  twenty  degrees  of  angle  of  attack. 
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